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Abstract 

In this paper, we present a new matrix approach for the analysis of subdivision schemes whose non- 
stationarity is due to linear dependency on parameters whose values vary in a compact set. Indeed, we 
show how to check the convergence in C £ (R S ) and determine the Holder regularity of such level and pa¬ 
rameter dependent schemes efficiently via the joint spectral radius approach. The efficiency of this method 
and the important role of the parameter dependency are demonstrated on several examples of subdivision 
schemes whose properties improve the properties of the corresponding stationary schemes. Moreover, we 
derive necessary criteria for a function to be generated by some level dependent scheme and, thus, expose 
the limitations of such schemes. 

Keywords: Level dependent (non-stationary) subdivision schemes, tension parameter, sum rules, Holder 
regularity, joint spectral radius. 


1. Introduction 


We analyze convergence and Holder regularity of multivariate level dependent (non-stationary) subdivi¬ 
sion schemes whose masks depend linearly on one or several parameters. For this type of schemes, which 
include well-known schemes with tension parameters 0, i El, I HR . ItR HI] , the theoretical results from [§] 
are applicable, but not always efficient. Indeed, if the level dependent parameters vary in some compact 
set, then the set of the so-called limit points (see [gj) of the corresponding sequence of non-stationary masks 
exists, but cannot be determined explicitly. This hinders the regularity analysis of such schemes. Thus, 
we present a different perspective on the results in Q and derive a new general method for convergence 
and regularity analysis of such level and parameter dependent schemes. The practical efficiency of this new 
method is illustrated on several examples. We also derive necessary criteria that allow us to describe the 
class of functions that can be generated by non-stationary subdivision schemes. Indeed, we show how to 
characterize such functions by the special property of the zeros of their Fourier transforms. 

Subdivision schemes are iterative algorithms for generating curves and surfaces from given control points 
of a mesh. They are easy to implement and intuitive in use. These and other nice mathematical properties 
of subdivision schemes motivate their popularity in applications, i.e. in modelling of freeform curves and 
surfaces, approximation and interpolation of functions, computer animation, signal and image processing 
etc. Non-stationary subdivision schemes extend the variety of different shapes generated by stationary 
subdivision. Indeed, the level dependency enables to generate new classes of functions such as exponential 
polynomials, exponential B-splines, etc. This gives a new impulse to development of subdivision schemes 
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and enlarges the scope of their applications, e.g. in biological imaging mm, geometric design [If], [f2| or 
isogeometric analysis SG3- 

The main challenges in the analysis of any subdivision scheme are its convergence (in various function 
spaces), the regularity of its limit functions and its generation and reproduction properties. The important 
role of the matrix approach for regularity analysis of stationary subdivision schemes is well-known. It allows 
to reduces the analysis to the computation or estimation of the joint spectral radius of the finite set of square 
matrices derived from the subdivision mask. Recent advances in the joint spectral radius computation Hm 
makes the matrix approach very precise and efficient. In the non-stationary setting, however, this approach 
has never been applied because of the several natural obstacles. First of all, the matrix products that 
emerge in the study of non-stationary schemes have a different form than those usually analyzed by the 
joint spectral radius techniques. Secondly, the masks of non-stationary schemes do not necessarily satisfy 
sum rules, which destroys the relation between the convergence of the scheme and spectral properties of its 
transition matrices. All those difficulties were put aside by the results in {sj, where the matrix approach 
was extended to general non-stationary setting. 

In this paper, in SectionjdJ we make the next step and consider level and parameter dependent subdivision 
schemes whose masks include tension parameters, used to control the properties of the subdivision limit. 
Mostly, the tension parameters are level dependent and influence the asymptotic behavior of the scheme. If 
this is the case, the scheme can be analyzed by Theorem 2], which states that the convergence and Holder 
regularity of any such non-stationary scheme depends on the joint spectral radius of the matrices generated 
by the so-called limit points of the sequence of level-dependent masks. In Theorem 3.5, we show that for 
the schemes with linear dependence on these parameters, the result of [§S, Theorem 2] can be simplified and 
be made more practical, see examples in Section [3. II In Section [|J we address the problem of reproduction 
property of subdivision schemes and of characterizing the functions that can be generated by non-stationary 
subdivision schemes. This question is crucial in many aspects. For instance, the reproduction of exponential 
pol ynomials is strictly connected to the approximation order of a subdivision scheme and to its regularity 
1191 ]. Essentially, the higher is the number of exponential polynomials that are being reproduced, the higher 
is the approximation order and the possible regularity of the corresponding scheme. 


2. Background 

Let M = ml £ Z sxs , \m\ > 2, be a dilation matrix and E = {0,..., |m| — 1} S be the set of the coset rep¬ 
resentatives of Z S /MZ S . We study subdivision schemes given by the sequence {5^(0, r £ N} of subdivision 
operators A a (r) : f(Z s ) —»• £(Z S ) that define the subdivision rules by 

(S a (r)c)(a) = a a-M/3 c (fi)’ a £Z S . 

The masks a (r ) = {a£ r \ a £ Z s }, r £ N, are sequences of real numbers ai r "* and are assumed to be all 
supported in {0,..., N} s , N £ N. For the given set 

OO 

K = Y J M~ l G , G = {-\m\,...,N + iy, (2.1) 

r=l 

the masks define the square matrices 

4 r) = (iMa+e-f,) „ r£N, S£E. (2.2) 

We assume that the level dependent symbols 

a {r \z) = £ a£V, = ^ Q1 • ■ • • • z?', z € (C\ {0}) s . 
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of the subdivision scheme 


c (r+1) = S a( r) c (r) = 5 a( r) S a( r- 1 ) ... S a( i) c (1) , r G N, 
satisfy sum rules of order 1+1, £ € No- For more details on sum rules see e.g Hi mg. 

Definition 2.1. Let l S No, r £ N. The symbol aS r >(z), z £ (C \ {0}) s , satisfies sum rules of order £ + 1 if 
a^(l,..., 1) = |m| s and max max \D' n a^ r \e)\ = 0, (2-3) 


where S = {e *T“T £ = (e *1™T £ 


hl<t e£H\{l} 


'i £s ), eel?} and D v = -... 


_gvs 

~dzT‘ 


The assumption that all symbols a^ r \z) satisfy sum rules of order £ + 1, guarantees that the matrices 
Ai r \ e e 1?, r e N, in El have common left-eigenvectors of the form 


(p(oO)aeif ’ P en <. 

where lb is the space of polynomials of degree less than or equal to £. Thus, the matrices A^\ e e E, 

r e N, possess a common linear subspace Vt C orthogonal to the span of the common left-eigenvectors 
(V) 

of A e , e G E, r G N. The spectral properties of the set 

T={4 r) k, eeE,re N} 

determine the regularity of the non-stationary scheme, see 0. 

Remark 2.2. In the univariate case, i.e. \m\ = \M\, the assumption that the symbols a^ r \z), r £ N, satisfy 
sum rules of order 1+1 implies that 

a ( r \z) = {l + z + ... + z^-y b^z a , z£C\ {0}, 

aez 


and 


4 r) k = (bflc+e-fi)' 




£ e E. 


(2.4) 


In the multivariate case, the explicit form of the matrices Ai r ^\v e , £ £ E, r £ N, depends on the choice of 
the basis ofVt, see e.g. ,3, Section 3.1] or jj]. 


Definition 2.3. A subdivision scheme {S a (r), r £ N} is C l - -convergent, if for any initial sequence c £ 
£oo(Z s ) there exists the limit function g c £ C £ (R S ) such that for any test function f £ C(R S ) 


lim 

k—> OO 


ffc(-) - S a(-)S a (r-D .. .S a d)c(a)f(M k • -a) 

Ol^lL s 


CL 


= 0. 


(2.5) 


For more details on test functions see 21]. Note that, it suffices to check (12.51) for only one test function /. 
Note also that, if all limits of a subdivision scheme belong to C^(R S ), then the scheme may not converge in 
C l , but only in C^R 8 ). 

In this paper, we also show how to estimate the Holder regularity of subdivision limits. 


Definition 2.4. The Holder regularity of the C°—convergent scheme {5 a (r>, r £ N} is a = £ + £, if £ is 
the largest integer such that g c £ C^(R S ) and £ is the supremum of v £ [0,1] such that 


max \D»g c (x) - D^g c (y)\ <\x- y \ v , 


x, y £ R s . 


We call a the Holder exponent o/{5 a (r), r £ N}. 
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The joint spectral radius of a set of square matrices was introduced in 41] and is independent of the choice 
of the matrix norm || • ||. 

Definition 2.5. The joint spectral radius (JSR) of a compact family A4 of square matrices is defined hy 

1/n 


p(M) := lim max 

n—>-oo Mi,..., M n £ M 


U M > 

j =i 


The link between the JSR and subdivision is well-known, see e.g. [bj, [H, 22, Sal SI- 


3. Parameter dependent subdivision schemes: matrix approach 

There are several examples of subdivision schemes that include a tension parameter. We call them param¬ 
eter dependent schemes. Often the tension parameter is level dependent and shows a certain asymptotic be¬ 
havior which implies the asymptotic behavior of the corresponding non-stationary scheme, i.e. lim a^ = a. 

r—too 

In this case, although the set {A^P, £ £ E, r £ N} is not compact, the convergence and regularity of the 
scheme {S^m, r £ N} can be analyzed via the joint spectral radius approach in Q. The results in [§J are 
still applicable even if the parameter values vary in some compact interval. Indeed, the existence of the limit 
points for the sequence {a^ r \ r £ N} of the subdivision masks is guaranteed, though these limit points are 
not always explicitly known. 

Definition 3.1. For the mask sequence {a^ r \ r £ N} we denote by A the set of its limit points, i.e. the 
set of masks a such that 


a £ A, if 3 {r n , n £ N} such that lim sS rn ' > = 


a. 


In this section, we show that the joint spectral radius approach can be effectively applied even if the 
limit points of {a^, r £ N} cannot be determined explicitly, but the masks a^ 7 ’) depend linearly on the 
parameter £ [wi,^], —oo < uj\ < u> 2 < oo. 

Well-known and celebrated examples of parameter dependent stationary subdivision schemes with linear 
dependence on the parameter are e.g. the univariate four point scheme [251 ] with the symbol 


a(z,uj) = ( ^ — htu(—z 2 + 1 + z 2 — z 4 ), co € 


°’T6 


z £ C\{0}, 


which is a parameter perturbation of the linear B-spline. Also the bivariate butterfly scheme 26j with the 
symbol 

a(zi,z 2 ,u>) = -(1 + Zi)(l + z 2 )(l + ziz 2 ) + u> c(z!,z 2 ), z lt z 2 £ C\{0}, 


with 


c(z i, z 2 ) = z 1 1 z 2 2 + z 2 z 1 1 + z 1 J, z 2 1 + zfz 2 1 - 2 zfz 2 - 2zfz 2 + zfz 2 + z\z 2 + z\z 2 

T zfz 2 — 2z 1 J- z^ — 2^ 2 — 2z 2 T zf T z 2 — 2^ 2 T z 2 (3.1) 

is a parameter perturbation of the linear three-directional box spline. Other examples of such parameter 
dependent schemes are those with symbols that are convex combinations 

ua(z!,z 2 ) + (1 - u) b(z lt z 2 ) = b(zi,z 2 ) +w(a(zi,z 2 ) - b(zi,z 2 )), w £ [0,1], zi,z 2 £ C\ {0}, 


of two (or more) symbols of stationary schemes, see e.g. [7j, [lj, |l6j, 29]. Known are also their non-stationary 
univariate counterparts with level dependent parameters (see i-itlSE , for example) 


^ + + ^ r) {~z~ 2 + 1 + z 2 - z 4 ), r £ N, lim w (r) = 


w £ 


2 
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and 


w (r) a(z) + (l-uj^)b{z), reN, w (r) e[0,l]. 

Note that the use of the level dependent parameters sometimes allows us to enhance the properties of the 
existing stationary schemes (e.g . with respect to their smoothness, size of their support or reproduction and 
generation properties 00 , 


laiia). 


In all schemes considered above, the subdivision rules depend either on the same, fixed, parameter u> = 
u £ \uj independent of r, or the parameters u>^ £ [wi,W 2 ] are chosen in a such a way that either 
lim = to £ [wi,W 2 ] or the corresponding non-stationary scheme is asymptotically equivalent to some 

r—>oo 

known stationary scheme. In this section, we provide a matrix method for analyzing regularity of more 
general subdivision schemes: we consider the level dependent masks a(u/ r )) = {a a (a/ r )), a £ Z s j, r £ N, 
and require that £ [uq, oq] without any further assumptions on the behavior of the sequence {o/ r ), r £ 
N}. We assume, however, that each of the masks depends linearly on the corresponding parameter u>( r ) m 
The level dependent masks {a(w ( ' r ^), r £ N} define the corresponding square matrices which we denote 

by 

(3.2) 


— (aMa+e-/3(w < ' r ' ) ) ) , £ £ E, 

\ / ql,p€lK 


and the level dependent symbols 


z(z,w w ) = ^2 a a (i 

a£Z‘ 


[U (r) )z° 


z£(c\m s 


The assumption that each mask a(u/ r )) depends linearly on cc- r ), leads to the following immediate, but 
crucial result. 

Proposition 3.2. Let £ £ No and —00 < uq < u >2 < 00 . If every symbol of the sequence {a(z, r £ N} 

depends linearly on the parameter uj^ £ [oq, oq] and satisfies sum rules of order £ + 1, then every matrix in 
T = {A e w (r) \v e , £ [uq, UJ 2 ] , £ £ E, r £ N} is a convex combination of the matrices withoj^ £ {aq,W 2 } 

= (1 -^ r) )A eiU11 \ Ve +t^A SiUl2 \ Ve , fW £ [0,1]. 

Proof. Let r £ N. We first write as a convex combination of uq and uq, i.e. 

w ( 1 ') = (1 — t^)uii + t^uj 2 with t ^ £ [0,1]. 

Note that all entries of the matrices A eul (r), £ £ E, are the coefficients of the corresponding mask a(cv^). 
Since the mask coefficients depend linearly on the parameter so do the matrices A e ul ( r ), and hence, 
the corresponding linear operators. Therefore, the restrictions of these operators to their common invariant 
subspace V) also depend linearly on this parameter. □ 

In the level independent case, i.e. = w for all r £ N, the use of the joint spectral radius approach 

for studying the convergence and regularity of the corresponding stationary subdivision schemes is well 
understood. To show how this approach can be applied in the our non-stationary setting, we need first to 
prove the following auxiliary result. 

Proposition 3.3. Let £ £ No and 

T = {A StU (.r) |vj, € [uq,uq], £ £ E, r £ N} (3.3) 

be the infinite family of square matrices. If the JSR of the family Tu, 1: u> 2 = {Ae, Ul \vi, Ar,oj a |vt> £ £ E} 
satisfies p(7Li,w 2 ) = 7 , then p(T) = 7 . 
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Pro of. F irst of all observe that p(T Ul W2 ) = 7 implies, for any 5 > 0, the existence of a (5-extremal norm (see 
e-g. HH), i.e. an operator norm || • ||,5 such that 

IIAr.ajilvf ||a < 7 + S, ||-^e,u 2 |vi II5 < 7 + (3-4) 

Then, by Proposition 13.21 estimates in (13.41) and subadditivity of matrix operator norms, we get 

IIAs.wMlvelk = II (1—|y e +^ r ^£,w 2 I Vi II (5 < (l-t (r) )||A £ , W i|vi||5+t (r) ||A £ , W2 |vi||5 = 7+^ f {r) € [0,1]- 

This, due to the arbitrary choice of S > 0, implies that p (T) = 7 , which concludes the proof. □ 

Remark 3.4. (i) Note that, if the family 7Li ,^ 2 ® s non-defective, i.e. there exists an extremal norm || • || 
such that max eE £; {II^E.ajilvilli ||A E|U , 2 \v t ||} = 7 > then T is also non-defective and all products of degree d 
of the associated product semigroup have maximal growth bounded by Note also that for any family of 
matrices B, B C T, it follows that p(B) < 7 . (ii) Moreover, if a family T is irreducible, i.e., its matrices 
do not have a common nontrivial subspace, then T is non-defective. Therefore, the case of non-defective 
families is quite general. 

We are now ready to formulate the main result of this section. 

Theorem 3.5. Let l £ No- Assume that every symbol of the sequence {a(z, u/ r )), r £ N} depends linearly 
on u/' r ) £ [oji, 012 } and satisfies sum rules of order £ + 1. Then the non-stationary scheme {S a ^ r )), r £ N} 
is C e -convergent, if the JSR of the family 7Li,w 2 = {A £!Ul \v e , A e ^ U2 \y l , s £ E} satisfies 

P(T Ul ,u 2 ) = 7 < H - *. (3-5) 


Moreover the Holder exponent of its limit functions is a > — log| m | 7 . 

Proof. Since the parameters {u/ r 4, r £ N} vary in the compact interval [wi, W 2 ], there exists a set of limits 
points (finite or infinite) for the sequence {a(a/ f 4), r £ N} of subdivision masks. Let us denote this set by A 
and the corresponding set of square matrices by 7a = {A e = {&Ma+e-p)a,pEK, £ £ E, a £ A}. Obviously, 
Ta C T with T as in (13.31) . Since by Proposition 13.31 and Remark [3Tl p(Ta) < 7 , the claim follows by [3, 
Corollary 1], □ 

Remark 3.6. (*) Note that, due to p(Ta) < 7 , Theorem \3 . 51 yields a smaller Holder exponent a than what 
could be obtained by Corollary 1]. For example, consider binary subdivision scheme with the symbols 


a(z,uj {r) ) = z r £ {1,..., L}, L £ N, 

fl( 2 ,w (r) ) = z " 1 ( 1 1 2) ~ + + z~ 1 +z- z 3 ), r>L + 1, 2 £C\{0}. 

z Id 

To apply Theorem 1,9.51 we can view the corresponding masks as being linearly dependent on parameters 
w ( ' r 4 £ [0, yg]. The corresponding family = {^ E ,o|vi! Ar.^lvi! £ £ {0,1}} consists of the four matrices 



/ 

— iO 

-2w + \ 

— LU 

0 

\ 


/ 

2 u> 

2 ui 

0 

0 

\ 


A),oj|vi = 


0 

2 co 

2cu 

0 

, A ltU1 \ Vl = 


—UI 

-2 w + i 

— Ul 

0 

(3.6) 


0 

—U! 

-2uj + \ 

—UI 



0 

2 

2co 

2uj 

0 



V 

0 

0 

2 u 

2 ui 

) 


V 

0 

— LU 

~2 co + \ 

—UI 

) 



for w £ {0, yg}. Due to 

max> { \\A Ei0 \ Vl \\oo, \\A e ^_ |v 7 ||°°} = {p(^,ok),P(-4 e ,ik)} = 


6 








we get p( 7q j_) = \ and, thus, the corresponding scheme is convergent and has the Holder exponent a > 1. 
On the other hand, the set A of limit points of the masks can be explicitly determined in this case and 
consists of the mask of the four point scheme. Thus, by fs, Corollary 1], the Holder exponent is actually 
a > 2 . 

(ii) The regularity estimate given in Theorem 1 3. 5\ can be improved, if the actual range of the parameters 
, r > L, for some L eN, is a subinterval of [lo\, 002 ], see section [Q1 

(in) Note that the result of Theorem \3.5\ is directly extendable to the case when the matrix family ff 
depends linearly on a convex polyhedral set fl = co{u>i,... ,u>l} of parameters u G fl C R p , r G N, such 
that 

L L 

with t^p G [0,1] and tP = 1. 
l=i l=i 

This is the case, for example, when we define the level and parameter dependent symbols 
a(z,u {r) ) =^2ujp aj (z), u > (r) = (ujP,...,lvP) t Gfl, r G N. 

1=1 


3.1. Examples 

In this section we present two univariate examples of level dependent parameter schemes, whose con¬ 
structions are based on the four point and six point Dubuc-Deslauriers schemes. In particular, in Example 
GUI the non-stationary scheme is constructed in such a way that the support of its limit function 


61 = lim S. 


a(uj ( r )) 


• • • 


1 , a = 0 , 

0 , otherwise 


a G T s , 


is smaller than the support of the four point scheme, but its regularity is comparable. In Example Id.81 every 
non-stationary mask is a convex combination of the four point and six point Dubuc-Deslauriers schemes. 
We show how the regularity of the corresponding non-stationary scheme depends on the range of the corre¬ 
sponding parameters {w^, r G N}. Both examples illustrate the importance of the dependency on several 
parameters r G N} instead of one u> G R. 

Example 3.7. We consider the univariate, binary scheme with the symbols 

a(zMp = rG {1,2}, 

a(z,^ r) ) = z- ^) +u) ( r )(-z~ 3 + z- 1 + z-z 3 ), r> 3, z€<C\{0}, 

where are chosen at random from the interval [^, ^]. The corresponding family 

%^ = {A efi \ Vl ,A e i_\ Vl , £ e {0,1}} 

consists of the same four matrices as in (1331) . And at the first glance the Holder exponent of this scheme is 
a > 1. On the other hand, we can view this scheme as the one with the corresponding matrix family 


Tjl jl 


applied to a different starting data. Then we get p(l~A. j _) = 3/8 and, by Theorem \3.5i the Holder exponent 

64 ’16 

is actually a > —log 2 f ~ 1.4150. 

The size of the support of (f\ can be determined using the technique from flS] and is given by 


_fc=0 k =0 


3 3 
2’2 
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with 


m 

m 


= — 1 , r(k) = 1 , k = 0 , 1 , 
= -3, r{k) =3, k > 2. 


Recall that the support of the basic limit function of the four point scheme is [—3,3]. 


Example 3.8. In this example we consider the univariate non-stationary scheme with symbols 
a{z,u} ( - r '>)=u;^a(z) + (l-^ r) )b{z), w (r) G [0,1], 2G<C\{0}, 

where 

, , Z~ 3 (Z + l) 4 / 2 a i'i 

«(*) =- yg( z -^ z + 1 ) 

is the symbol of the four point scheme and 


b(z) = — — ^ ^ (3 z i — 18z 3 + 38z 2 — 18z + 3) 
256 v ' 


is the symbol of C 2 —convergent quintic Dubuc-Deslauriers subdivision scheme fUQ]. By fH/, the Holder 

exponent of the Sh is a 


2.8301. 

To determine the regularity of this level and 

parameter dependent scheme 

we consider the matrix set 











75,1 = {71 

:,o|v 2 !^e,l|v2> £ e (0; 1}} 




with the four matrices 

/ 

3 — 3w 

0 

0 

0 

0 

0 

0 

\ 




-7-9 u 

—9 + 9 lo 

3-3w 

0 

0 

0 

0 



. . 1 

A,, w |y 2 - 256 


45 + 3w 

45 + 3 uj 

-7-9 w 

—9 + 9 uj 

3 — 3w 

0 

0 




—9 + 9w 

-7-9 u 

45 + 3w 

45 + 3w 

-7-9 u 

—9 + 9w 

3 — 3w 
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0 

3 — 3w 

—9 + 9 uj 

-7- 9w 

45 + 3w 

45 + 3 w 

-7-9 u 





0 

0 

0 

3 — 3w 

—9 + 9uj 

-7-9 u 

45 + 3 w 




V 

0 

0 

0 

0 

0 

3 — 3w 

—9 + 9w 

) 



/ 

—9 + 9 to 

3 — 3 uj 

0 

0 

0 

0 

0 

\ 




45 + 3w 

-7-9 u 

—9 + 9 u> 

3- 3w0 

0 

0 




Al ’ uWa = 2^6 


-7-9 to 

45 + 3w 

45 + 3w 

-7-9 uj 

—9 + 9 lo 

3 — 3w 

0 




3 — 3w 

-9 + 9w 

-7-9 w 

45 + 3 w 

45 + 3 uj 

-7-9 u 

—9 + 9w 




0 

0 

3 — 3cv 

—9 + 9w 

-7-9 u 

45 + 3w 

45 + 3 ui 





0 

0 

0 

0 

3 — 3w 

—9 + 9 w 

-7-9 w 




V 

0 

0 

0 

0 

0 

0 

3 — 3w 

) 



for u) G {0,1}. In this case, the regularity of the non-stationary scheme {S^^Mp r G N} coincides with the 
regularity of the four point scheme. For w G {a, 1}, a > 0, the scheme {£„(„(»•) p r G N} is C 2 —convergent. 
And, for uj G {0, a} , a < 1, extensive numerical experiments show that the JSR of the family 7o,a is 
determined by the the subfamily {A e>a \y 2 , £ G {0,1}}. For example, for a = we obtain p(T 0 i) ~ 0.2078 
and, thus, the corresponding Holder exponent is a > 2.2662. 


4. Limitations of generation properties of non-stationary schemes 


It is known that certain level dependent (non-stationary) subdivision schemes are capable of general 


ing/reproducing certain spaces of exponential polynomials, see e.g. [9|, 18]. In this section, we are interested 


in answering the question: How big is the class of functions that can be generated/reproduced by such 
schemes? 









More precisely, we show that, already in the univariate setting, the zero sets of the Fourier transforms of 
the limit functions 

= <(«> = {;; otherwise ’ “ £ Z ‘' 

of such schemes are unions of the sets 

IV = (w e C : a (r )(e- i2,rM “ r “) = 0}, r > k, 

and that the sets T r are such that r r + M rr L = P r . Thus, some elementary functions cannot be generated 
by non-stationary schemes, see Example 14.41 Also the requirement that 

cj) k (oj) = f (f k (x)e~ l2nxui dx, w G C, k e N, 

Jr 

is an entire function, limits the properties of the functions that can be generated by non-stationary subdi¬ 
vision schemes. 

Proposition 4.1. Let {<f k , k € N} be continuous functions of compact support satisfying 

4>k(x) = E a( fc )(a)0fc+i (Mx-a), ke N, 

ol&L 


Then 


{w e C : Sk{uj) = 0} = 1J T r , 

r>k 


such that the sets T r satisfy 


r r + M r z = r r . 


Proof. Let k € N. By Paley-Wiener theorem, the Fourier transform (j> k defined on M has an analytic 
extension 

4>k(v) = [ 4> k {x)e~ l2 ' KXU dx , weC, 

Jr 

to the whole complex plane C and i k is an entire function. By Weierstrass theorem j20 ], every entire function 
can be represented by a product involving its zeroes. Define the sets 

r r := {w e C : r «) = 0 }, r £ N. 

Let z r p ,..., z r ,N be the zeros of the polynomials a^ r \e~ l2 ^ M “), counting their multiplicities. Then 

N 

r r = iM r (J Ln (z rt e), 


where, by the properties of the complex logarithm, each of the sets iM r Ln(z ri i) consists of sequences of 
complex numbers and is M r — periodic. Thus, each of the sets T r satisfy 


T r + M r Z = T r , re N. 

The definition of <j> k as an infinite product of the trigonometric polynomials a ( - r \e~ l2nM “), r > k, yields 
the claim. 

□ 


The following examples illustrate the result of Proposition 14.11 
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Example 4.2. The basic limit function of the simplest stationary scheme is given by cj>\ = X[o.i)- Its Fourier 
transform is 

„ 1 — e ~i2Trui 

(j>i (w) = ---, and {u> £ C : cj>i(uj) = 0} = Z \ {0}. 

iznuj 

The mask symbol a{z) = 1 + z has a single zero at z = — 1, i.e. e -12 ’ 1 ’ 2 “ = —1 for uj = 2 r {i + k k £ Z}, 
r £ No- In other words, Ti = {1 + 2fc : fcgZ} and r r = 2r r _i for r > 2. Therefore, 

{w £ C : <fi (cv) = 0 } = (J T r . 

r£N 


Example 4.3. The first basic limit function of the simplest non-stationary scheme is given by 4>\{x) = 
X[o,i)( x ) eXx > ^ G C. Its Fourier transform is 

„ p -i2nu>+\ _ i » i \ 

= -—- —OJ € C, and {w £ C : <fi(co) = 0} = -— + Z \ {0}. 

—IZ7TU) + A Z7T 

The mask symbol a( r \z) = 1 + e A2 z has a single zero at z = —e _A2 , i.e. e ~ l2w2 “ = —e _A2 for 
w = — + 2 r {i + k : k £ Z}, r £ N. Note that Ti = — |^ + {l + 2fc : fc£Z} and 

(J 2 r {\ + k : fc £ Z} = Z \ {0}. 

reN 


Therefore, 

{w £ C : </>i (w) = 0 } = [J T r . 

reN 

In the next example we identify a compactly supported function that cannot be generated by any non¬ 
stationary subdivision scheme. 

Example 4.4. Let us consider the compactly supported function 

2 

/(ff) =%[-!,!]( x) - =, z£R. 

V 1 — £ 

It cannot be a limit of any non-stationary subdivision scheme. Indeed, its Fourier transform 

J 0 (u) = [ f{x)e- ixu dx, oj £ C, (4.1) 

Jr 

is the Bessel function Jo of the first kind, which is entire, but has only positive zeros. The lower bound for 
its zeros jo jS , s £ N, is given by jo tS > \J {s — see [36]. Thus, Proposition \f.l\ implies the claim. 
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